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Abstract 

Plastidic uracil salvage Is essential for plant growth and development. So far, PLUTO, the plastidic nucleobase transporter 
from Arabidopsis tlialiana is the only known uracil importer at the inner plastidic membrane which represents the 
permeability barrier of this organelle. We present the first homology model of PLUTO, the sole plant NCSl member from 
Arabidopsis based on the crystal structure of the benzyl hydantoin transporter MHPl from IVIicrobacterium liquefaciens and 
validated by molecular dynamics simulations. Polar side chains of residues Glu-227 and backbones of Val-145, Gly-147 and 
Thr-425 are proposed to form the binding site for the three PLUTO substrates uracil, adenine and guanine. Mutational 
analysis and competition studies identified Glu-227 as an important residue for uracil and to a lesser extent for guanine 
transport. A differential response in substrate transport was apparent with PLUTO double mutants E2270 G1470 and E227Q 
T425A, both of which most strongly affected adenine transport, and in V145A G1470, which markedly affected guanine 
transport. These differences could be explained by docking studies, showing that uracil and guanine exhibit a similar 
binding mode whereas adenine binds deep into the catalytic pocket of PLUTO. Furthermore, competition studies confirmed 
these results. The present study defines the molecular determinants for PLUTO substrate binding and demonstrates key 
differences in structure-function relations between PLUTO and other NCSl family members. 
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Introduction 

Salvage of nucleosides and nucleobases represents an alternative 
to the de novo synthesis of nucleotides. Hereby, energy and also high 
valuable constituents of purine and pyrimidines for plant 
development namely phosphate and nitrogen are preserved. It 
was shown that salvage of uracil and of uridine in plastids is 
essential for growth and development of Arabidopsis seedlings 
[1,2]. Furthermore, stimulation of the pyrimidine salvage pathway 
in potato leads to an increase in biosynthetic performance and 
starch contents in tubers [3]. 

Our recent results indicate that pyrimidine de novo synthesis is 
finalized in the cytosol and consequently pyrimidines must be 
initially imported into plastids to serve for the synthesis of nucleic 
acids in this compartment by salvage reactions [4]. Whereas no 
nucleoside transport system has been described for plastids so far, 
a corresponding nucleobase transporter was identified. PLUTO, 
the plastidic nucleobase transporter is the sole member of the 
NCSl (Nucleobase:Cation Symporterl) protein family in Arabi- 



dopsis [4,5]. Other plant species like rice {Oryza sativd), maize [Z^a 
mays), wine {Vitis viniferd) or poplar (Populus trichocarpd) similarly only 
harbor one PLUTO homolog. However, in the grass species 
Bravhjpodium distachyon two PLUTO homologs are present [6]. 
Recently, a PLUTO homolog from the single cell algae 
Chlamydomonas reinhardtii was identified, mediating transport of 
adenine, guanine, uracil and aUantoin. This shows that solute 
specificity of plant NCSl members arose early in plant evolution 
[7]. 

PLUTO plant homologs show an N-terminal extension counted 
from the beginning of the first transmembrane span of about 80 
amino acids. From the species mentioned above, only P. trichocarpa 
exhibits a shorter N-terminus (38 amino acids), whereas 
Arabidopsis PLUTO exhibits a much longer N-terminal extension 
of 134 amino acids [6]. Such sequence extensions are typical for 
eukaryotic transport proteins [8] and in fact prokaryotic NCSl 
family members like MHPl, HYUP and PUCI do not show such 
extensions [4] . Heterologous expression of the full length PLUTO 
Arabidopsis gene in Escherichia coli uroA insertion mutants 
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demonstrated PLUTO facilitates high affinity nucleobase trans- 
port. The apparent affinities of uracil, adenine and guanine were 
determined with 16.4, 0.4, and 6.3 |iM, respectively [4]. While fuU 
length PLUTO can transport nucleobases, protein synthesized 
from a second putative start codon at position +97, almost lacking 
the complete N-terminal extension was not functional in the E. coli 
system and showed no targeting to chloroplasts in transient 
expression studies [4]. However, in the EMBL database PLUTO 
protein is annotated as N-terminal truncated ORF (Q93726) and 
was shown to function better in a yeast expression system, 
compared to the longer ORF [5]. At present it is unknown 
whether both protein variants exist in Arabidopsis. 

Plastids are plant specific organelles fulfilling a multitude of 
biochemical reactions such as photosynthesis, nitrogen and sulfur 
assimilation, starch, fatty acid and amino acid synthesis. For the 
necessary exchange of pathway intermediates transport proteins 
are located in the inner envelope membrane of plastids, which 
serves as a permeability barrier in these organelles. Amongst these, 
transporters for triose-phosphates, sugars and phosphate have 
been characterized (for review see [9]). In addition to PLUTO 
further transporters catalyzing the exchange of intermediates of 
nucleotide metabolism exist in plastids. These facilitate the 
exchange of ATP against ADP (NTT, [10,11]) and the net 
exchange of ATP, ADP and AMP (BTl, [12]). 

NCSl proteins like PLUTO are predicted to contain 12 
transmembrane spans and are also known as purine-related 
transporters. Being symporters they generally transport purines 
using a co-transport with protons. Apart from plants they can be 
found in archaea, bacteria, yeast and fungi [13] (www.tcdb.org). 
Well-characterized NCSl proteins with marked homology to 
PLUTO are the uracil transporter FUR4 and the uridine 
transporter FUIl from Saccharomyces cerevisiae [14,15], the allantoin 
transporter DAL4 from Saccharomyces cerevisiae [16,17] and FCYB 
from Aspergillus nidulans. The latter was characterized as purine- 
cytosine transporter and analyzed with respect to structure- 
function relationships recently [8]. 

The first structural information about the NCSl transporters 
came through MHPl, the benzyl-hydantoin transporter from 
Microhacterium liquefaciens. The MHPl crystal structure revealed a 
monomer organized into an alpha-helix bundle with 12 trans- 
membrane (TM) helices that form a deep internal cavity open on 
the periplasmic side with the substrate binding site located at the 
center of the protein cavity. The 12 TMs are arranged in two 
repeating halves (TMs 1 to 5 and TMs 6 to 10), connected by an 
18-residue loop and followed by two additional TMs. Qiiite 
remarkably, MHPl structure has been solved at different stages 
along the transport cycle, inward facing, occluded and outward 
facing states, which provides the structural basis for substrate 
transport [18]. By the analysis of these three MHPl conformations 
coupled with molecular dynamics simulations a transport mech- 
anism for MHPl has been proposed [19]. The conformational 
transition from the outward state to the inward state results from 
rigid body motion of the TM hehces 3, 4, 8 and 9 relative to the 
rest of the protein. This provides the molecular basis for the 
alternative access mechanism of substrate transport, which could 
be possibly apphcable to all the members of NCSl family. This 
"rocker switch" mechanism ensures "uphill" transport of substrate 
against its concentration gradient when it is coupled to energet- 
ically favorable "downhill" transport of the driving ion Na"''. 

The MHPl complex with benzyl hydantoin along with other 
conformations provides insight into the substrate transport 
mechanism. However, details about the substrate binding sites of 
NCSl family members for other known substrates including 
purines, cytosine and uracil are lacking. Despite a common 



predicted structural fold with the NCSl family, these transporters 
share divergent amino acid sequences which may lead to 
variations within the substrate binding sites and transport of a 
wide variety of substrates. In this study, we aim to address this 
issue by identifying critical residues involved in the binding of the 
physiological substrates namely adenine, guanine and uracil within 
the sole plant NCSl family member PLUTO. We provide an 
experimentally validated three-dimensional model of the binding 
site in PLUTO using mutational analysis of PLUTO paired with 
structural-activity data for the substrates. Validation of binding 
models was performed based on the mutational dataset of PLUTO 
point mutants describing the contribution of residues in the 
putative substrate binding pocket. Together with results from 
substrate competition studies mechanistic insights into PLUTO 
transport function were obtained. 

Materials and Methods 

Heterologous expression of PLUTO and site directed 
mutants in £ coli JD23420 

For the uptake experiments, PLUTO was heterologously 
expressed in an Escherichia coli transposon insertion strain lacking 
the endogenous uracil permease uraA (JD23420), obtained from 
the National Institute of Genetics (Shizuoka, Japan; http://www. 
shigen.nig.ac.jp/ecoh/strain/top/top.jsp). PLUTO and site di- 
rected mutants were cloned in the bacterial expression vector 
pTAC-MAT-Tag-2 (Sigma Aldrich, Deisenhofen, Germany) as 
previously described [4] . Aft<;r the transformation of the cells with 
the PLUTO expression plasmids, the cells were grown at 37°C and 
1 70 rpm in YT medium containing ampicillin and kanamycin. 
The expression of PLUTO and site directed mutants was initiated 
by the addition of isopropyl P-D-thiogalactopyranosid (IPTG) with 
a final concentration of 0.02% at an ODgoo of 0.5. The cells were 
grown for 2 h after induction, collected by centrifugation for 
5 min at 4500 g and 4°C and resuspended with potassium 
phosphate buffer (50 mM, pH 7.0) to an ODgoo of 5. The cells 
were stored on ice until use. 

Site directed mutagenesis 

Site directed mutagenesis was achieved by using the Quick- 
Change II-E Site directed mutagenesis kit (Stratagene) according 
to the manufacturer's instructions. Primers used are listed in Table 
SI. 

Nucleobase uptake experiments with PLUTO and site 
directed mutants 

For the uptake experiments with radiolabelled nucleobases, 
100 (jl IPTG-induced cells and uninduced cells used as negative 
controls were added to the same volume of preheated potassium 
phosphate buffer (50 mM, pH 7.0, 37°C) containing radiolabelled 
nucleobases (['*C]-uracil, -guanine and -adenine, 1012 Bq 
mol~'). The uptake of nucleobases with concentrations of 
20 [iM was performed at 37°C and permanent shaking at 200 
rpm. For the analysis of PLUTO mutants the cells were incubated 
for 5 min (2 min for competition studies) and the uptake was 
terminated by vacuum filtration through a membrane filter 
(0.45 |lm pore size; Whatman). The' cells were washed three times 
with 1 ml ice cold potassium phosphate buffer (50 mM, pH 7.0) to 
remove unimported nucleobases and the radioactivity was 
quantified in a Packard Tricarb 2500 scintillation counter after 
the addition of 4 mL scintillation cocktail (Roth). For the 
competition studies, all competitors were added with ten-fold 
excess (200 (tM) to the uptake medium containing 20 |.iM 
radiolabelled nucleobases. To show that PLUTO and PLUTO 
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mutants have no marked differences in protein amount, 
membranes from sonified E. coli cells were collected by sequential 
centrifugation, spotted on hybridization membranes and probed 
with a MAT-Tag specific antibody (Figure SI). 

Construction of the Sequence Alignment 

Multiple alignments of protein sequences were performed using 
ClustalW [20] and visualized with Genedoc. Sequences used are 
as follows: Arabidop.m thaliana PLUTO (At5g03555, AED90625), 
Microbacterium liquefaciens MHPl (2JLN_A) Aspergillus nidulans FCYB 
(Gl: 169798762), PLUTO homologs from Oryza saliva (Os02g4 
4680), Zea mays (Zm362848), Vitis vinifera (GSVFVTOIOSS 
705001), Populus tnchocarpa (Pt0006S121 10), and Brachypodium 
distachyon (Bradi3g51350). 

Homology modeling 

The three-dimensional model of the PLUTO full length 
sequence was built using crystal structures of MHPl in the 
outward-facing open conformation (PDB ID 2JLN) as a template. 
Structural modeling was performed using the 1-TASSER server 
based on ah inrfzo/ threading method [21]. The server predicted 5 
models of PLUTO and the best model was selected based on C- 
score (0.59). C-score is a confidence score for estimating the 
quality of predicted models by I-TASSER. It is calculated based 
on the significance of threading template alignments and the 
convergence parameters of the structure assembly simulations. 
Quality control analysis by the PROCHECK program identified 
90% residues of the model in the most favorable regions of the 
Ramachandran plot [22]. Only five glycine and nine non-glycine, 
non-proline residues (1% and 1.8% of the primary structure) were 
in the forbiddcm region of the Ramachandran plot. For the 
evolutionary conservation profile of the PLUTO core protein 
(TMI-TMIO) the conservation score was calculated based on 
ClustalW matrix by Conserv Server [23]. 

Molecular dynamics simulations 

The initial homology model was subjected to protein prepara- 
tion tool implemented in Schrodinger suite 2012. The N-termincd 
part of PLUTO which does not align to MHPl was removed. In 
summary the hydrogen atoms were added and the protonation 
states of histidine residues were optimized based on pKa 
calculation. A restrained minimization was performed using 
Impact MM engine implemented in Schrodinger suite 2012. 
The molecular dynamics (MD) simulation was performed using 
Desmond v.3 software. The protein structure was embedded into a 
l-palmitoyl-2-oleylphosphatidylcholine (POPC) lipid bilayer. The 
length and width of lipid box was 105 A x 1 1 7 A. The system was 
solvated using TIP4P explicit water and neutrahzed by adding 
150 mM sodium ions. After membrane insertion and solvation of 
the system, the lipids were re-equilibrated for 10 ns. The MHPl 
template available in OPM database was used to set orientation of 
protein in the membrane [24]. At last, a total of 20 ns MD 
simulations were performed on th(" system under a constant 
temperature of 300 K and a constant pressure of 1 bar. 

Molecular Docking 

The protein structure obtained after homology modeling and 
MD simulations were employed in docking calculation performed 
with Glide 5.0 (Schrodinger, LLC). AU three substrates were 
minimized using the OPLS-2005 force field by Macro Model 
software [25] . Details of the algorithm are available in the GLIDE 
documentation. Briefly, GLIDE proprietary conformational 
expansion and exhaustive search of the binding site produce a 



multitude of ligand poses, which undergo an initial refinement, 
energy minimization on a precomputed grid, and a final scoring 
and ranking. GLIDE uses proprietary scoring functions that are 
variations of the Chem Score 43 empirical stxjring function and 
the OPLS-2005 force field to compute van der Waals and 
electrostatic grids for the receptor. The final ligand binding poses 
were ranked according to a computed Emodel score that 
encompasses the grid score, the proprietary Glide score, and the 
internal energy strain. The substrates were docked into the 
receptor site using GLIDE SP (standard precision) mode. A total of 
50 docked conformers for each substrate were obtained. The best 
pose for each ligand was selected using total energy as a primary 
criterion and the experimental knowledge as a secondary criterion. 

Results 

Identification of critical residues for substrate binding 
through sequence homology 

PLUTO (At5g03555), identified as plastidic nucleobase trans- 
porter in Arabidopsis [4], belongs to the Nucleobase:Cation 
Symporterl (NCSl) family of transport proteins. Further NCSl 
members are FUR4, FUI 1 , and DAL4 from Saccharomyces cerevisiae, 
HYUP from Arlhrobacler aurescens, PUCI from Bacillus subtilis, and 
MHPl from Microbacterium liquefaciens. Among these sequences 34 
amino acid residues are highly conserved including six residues 
known to be involved in substrate binding based on the available 
crystal structure of MHPl [26]. From sequence alignment with 
MHPl, FCYB and six further sequences from plants (Figure lA), 
we observe that four of these residues, Trp-223, Trp-342, Asn-426 
and Asn-430, relative to PLUTO are conserved in all sequences, 
whereas two, Gly-147 and Glu-227 difier. Gly-147 of PLUTO is 
equivalent to Gln-42 in MHPl and Glu-227 correspondingly 
aligns to Gin- 121 [4] . All six named positions are conser\-cd in all 
plant species shown in the alignment except for Vitis vinifera 
possessing a G147T exchange (Figure lA). Based on the 
information from the ahgnments we aimed to gain insight into 
the structure function relationship of nucleobase binding and 
transport mediated by PLUTO. To do so, both non conserved 
amino acid residues in PLUTO, Gly-147 and Glu-227 which are 
involved in substrate binding in MHPl were mutated and the 
protein was subserjuently tested for transport activity after 
expression in E. coli cells as previously described [4]. Heterolo- 
gously expressed PLUTO carries a C-terminal MAT-Tag (Sigma 
Aldrich) used to identify the protein in western blots (Figure SI). 
The tag does not negatively influence transport activity of 
PLUTO. In contrast, the background of uninduced controls was 
lower when PLUTO was synthesized with tag (Figure S2). 

Functional studies with PLUTO mutants 

In the mutagenesis approach, first, Glu-227 was exchanged to 
an aspartate residue also carrying a negative charge. E227D 
showed markedly reduced transport of uracil and guanine (Table 1, 
Figure 2). Whereas wild type PLUTO protein catalyzed import of 
3.7±0.28 nmol ['*C]-uracil mg"' protein (20 llM) into induced 
E. coli cells, uninduced cells, used as negative controls, only 
imported 0.5±0.11 nmol [' 'C] -uracil mg ' protein within 5 min 
incubation leading to a net import of 3.2 nmol ['*C] -uracil mg ' 
protein (Figure 2A). In E227D expressing cells uracil transport was 
reduced to 0.53 nmol mg ' protein. Guanine transport (20 pM) 
was reduced to 0.86 nmol mg ' protein relative to controls which 
accounted for 3.4 nmol mg~' protein (Figure 2B). In contrast, 
adenine transport (20 pM) was only slightly reduced to 2.1 nmol 
mg ' protein representing 64% of the transport rates of controls 
(Figure 2C). The calculated values in this mutant were 5.18, 
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Figure 1. Alignment of PLUTO witKi otiier NCSI-type protein sequences from different organisms and PLUTO liomoiogy modei. (A) 

NCSl proteins from Arabidopsis thaliana PLUTO (AED90625), Microbacterium liquefaciens MHP^ (2JLN_A) Aspergillus nidulans FCYB (Gl: 169798762), 
and PLUTO homologs from Oryza sativa (Os02g44680), Zea mays (Zm362848), Vitis vinlfera (GSVIVTOl 033705001), Populus trichocarpa 
(Pt0006S12110), Brachypodium distachyon (Bradl3g51350) were aligned with ClustalW (www.ebi.ac.uk). The residues are shown in white, gray or 
black color according to their conservation mode (black for highly conserved residues). Residues marked with asterisks were mutated in the course of 
this work. Residues marked with red boxes are directly involved in PLUTO substrate binding. Blue boxes indicate two Trp residues which are highly 
conserved among NCSl proteins and probably exhibit a function in stabilizing the protein-substrate complex with weak pi-stacking interactions. (B) 
PLUTO homology model with marked TIVIs. A three-dimensional model of PLUTO was built using l-TASSER server based on structural information of 
MHPI. The TIMs are marked in colors according to the alignment in Figure 1A. The N-terminus of PLUTO was removed for better visualization. 
doi:l 0.1 371/journal.pone.0091 343.g001 



3.70 and 0.97 for uracil, guanine and adenine, respectively 
(Table 2). In the next step negatively charged glutamate at position 
227 was mutated to an uncharged glutamine (E227Q), a smaller 
and uncharged alanine (E227A), or to a positively charged lysine, 
E227K. In any of these cases uracil transport activity was 
completely extinguished (Figure 2A) and guanine transport was 
reduced to 21%, 4% and 11%, respectively (Figure 2B). However, 
adenine transport even increased in E227Q^by 21% (Figure 2C) 
whereas the affinity remained nearly unchanged (Table 2). E227A 
and E227K mutations led to residual adenine transport of 26% 
and 7% compared to controls, respectively (Figure 2C). Wild t}pe 
and mutated PLUTO expressing cells showed no marked 
differences in protein amount. For this, membranes from sonified 
E. coli cells were collected by sequential centrifugation, spotted on 
membranes and probed with a MAT-Tag specific antibody 
(Figure SI). No marked differences in protein amount were 
observed. 

Mutation of the second site of interest based on the sequence 
comparison between PLUTO and MHPI, Gly-147 was also 
carried out. Adenine transport was impaired by 45% when Gly- 
147 was mutated to glutamine G147Qj whereas uracil and 
guanine transport were not significandy altered (Figure 2 A-C). 
The affinity towards uracil in this mutant slightly decreased 
whereas affinities for adenine and guanine remained nearly 
unchanged (Table 2). However, in the double mutant E227Q, 
G147Q, which resembles the MHPI structure with respect to the 
identified substrate binding site residues, transport of eiU the 
identified substrates uracil, guanine and adenine was reduced to 
16% or below (Figure 2 A-C, Table 1). In contrast to the single 
mutation E227Q_also adenine was no longer transported although 
the protein amount in E. coli membranes was not changed 
compared to controls and this holds true for all other mutants 
described in this paragraph (Figure SI). 

PLUTO structural model 

To obtain further insight into the PLUTO structure-function 
relationship, wc took advantage of the close sequence homology 
between PLUTO and MHPI and developed a homology model 
based on the outward facing conformation of MHPI [18] (Fig IB, 
3A). The overall protein sequence homology between MHPI and 
PLUTO (26% sequence identity and 48% sequence similarity) is 
adequate for sustaining a theoretical model [21]. The model 
building was performed using I-TASSER and 5 top-scoring 
homology models of PLUTO were built. The stereochemical 
quality of the models was further analyzed using PROCHECK 
[22]. In the best scoring model, unfavorable protein conformations 
were limited to regions for which no template structures were 
available. In the next step, we performed MD simulations on the 
PLUTO model for 20-ns to measure stability of the structure. 
Therefore, the protein model was embedded in POPC lipid 
bilayer and explicit water was utilized for solvation of the 
membrane. During the MD simulations RMSD of TM helices 
Coi atoms from starting coordinates were monitored. Throughout 
the simulations RMSD of Ca atoms were steady and did not 



exceed more than 3 A, suggesting the stability of the homology 
model (Figure S3). 

The structural model of PLUTO (for PDB file see File SI) is 
composed of 12 TM helices, in which 10 TMs form a compact 
core foUowed by two additional TMs (Fig IB, 3A). The PLUTO 
structure contains the characteristic '5-H5' inverted structural 
symmetry motif defined by TMI-TMIO. Within the "inverted 
repeat" motif, two sets of five TMs are oppositely oriented with 
respect to the each other around an axis in the center of the 
membrane and parallel to the plane. For convenience, the 
structure can be divided into two topologically distinct sub- 
domains. The first is a four helix bundle comprising TMs 1 and 2 
and their pseudo two-fold equivalents TMs 6 and 7 are proposed 
to form the substrate binding site (Figure 4A). The second is 
another motif of four helices formed from TMs 3 and 4 and their 
pseudo two-fold equivalents, TMs 8 and 9 and forms a possible 
ion-binding site [26] . Conserved residues are located in the core of 
the structure (Figure 3B). The substrate binding site is located in 
the center of the helical bundle (Figure 3A, 4A). 

Substrate docking 

All the physiological substrates of PLUTO namely adenine, 
guanine and uracil were docked into the substrate binding pocket 
of the protein using the standard precision (SP) protocol. All of 
these substrates bind into a site located in the middle of the 
membrane-spanning regions, where they form direct contact with 
residues from TMl, TM3, TM6 and TMS (Figure 3A, 4A). 
Residues in these regions have also been found important for 
ligand binding in MHPI [26] and FCYB [8], which suggest that 
the fold and function of this protein region is highly conserved 
among the NCSl members. This is also supported by an 
evolutionary conservation profile as shown (Figure 3B). Guanine 
interacts with PLUTO through a bidentate H bond between 
purine N9 and residue Glu-227 and two H bonds formed between 
purine C2-NH2 and backbone carboxyl groups of residue Val-145 
and Gly-147 (Figure 4B), uracil adopts a binding mode in which 
pyrimidine N3 of uracil forms a H bond with residue Glu-227 
(Figure 4C). While both guanine and uracil interact with Glu-227 
through H bonds, in contrast adenine binds deep into the catalytic 
pocket of PLUTO located in a 5 A distance from the side chain of 
Glu-227 (Figure 4D). Adenine interaction is stabilized by three H 
bonds, one between purine N7 and Val-145 and a bidentate H 
bond between purine C6-NH2 and carboxyl group of Thr-425. In 
addition adenine exhibited weak /if-stacking interaction with 
indole rings of Trp-342. 

Validation of docking model by further functional studies 

Using docking studies as an initial guide we identified further 
critical rc-sidues that are located within 5 A of the putative binding 
site in addition to Gly-147 and Glu-227. We substituted each 
position with one or more residues having different physiochem- 
ical properties of the wild-type residue. Residues in the proximity 
of the proposed binding pocket that were mutated in the course of 
this work and shown to affect transport properties are listed in 
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Table 1. Nucleobase transport activity of different PLUTO 
mutants. 



Mutant 



Uracil transport Guanine transport Adenine transport 



E227Q 
E227A 
E227K 
E227D 
G147Q 
G147K 

E227Q G147Q 
E227Q T425A 
E227Q V145A 
V145A 

V145A G147Q 

T425A 

W342A 

W223A 

L144A 

F341A 

N426A 

N430A 

I226A 



+ 
++ 



The transport of uracil, guanine and adenine, 20 each, was measured after 
heterologous PLUTO expression in E. coli cells lacking the endogenous uracil 
importer uraA. Transport was corrected for uptal<e into uninduced cells. 
Transport function above 80% of control is indicated with "+++ ", between 80% 
and 50% with "++", between 50%and 20% with "+" and below 20% with 
doi:l 0.1 371/journal.pone.0091 343.t001 



Table 1 and their position is sliown in Figure I A. The transport of 
all tlie PLUTO substrates in all the single mutants, VI 45 A, 
T425A, W342A, and W223A, L144A, F341A, N426A, N430A 
and I226A was not markedly altered (Table 1). However, 
interesting observations for the double mutants of critical residues 
were made. The single mutant T425A had no effect on the 
transport but when a double mutant E227Q,T425A was tested, it 
reduced guanine and adenine transport to 25% and 32%, 
respectively (Figure 2B, C). In the docking pose both guanine 
and adenine appear to form a H bond with the backbone carboxyl 
group of Thr-425, however uracil does not show any interaction 
with Thr-425. Mutation of Thr-425 with any other amino acid 
residue will change the side chain only but the backbone carboxyl 
group win remain intact, which probably explains why the single 
mutant with Thr-425 had no effect on the transport. However, 
when the double mutant E227Q^T425A was studied the effect of 
Thr-425 binding was more visible, especially on adenine transport, 
which was shown to be unaffected by the single mutants E227Q, 
and T425A. A similar trend, although to a lesser extent was 
observed with double mutant V145A G147Qj which reduced 
guanine transport to 12% relative to controls. Interestingly, uracil 
transport in this mutant was even increased by 50% (Figure 2A). 
The mutational analysis of the putative substrate binding site 
revealed that uracil transport is affected by Glu-227, guanine 
transport by Glu-227, Thr-425, Val-145 and Gly-147 whereas 
adenine transport reUes on Thr-425 and Gly-147. 
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Figure 2. Functional analysis of PLUTO mutants. Radiolabeled 
nucleobase uptake studies with all substrates of PLUTO, namely uracil, 
guanine and adenine were performed after heterologous PLUTO 
expression in £ coli cells lacking the endogenous uracil transporter 
uraA. Uptake was measured for 5 min as previously described in [3] 
with nucleobase concentrations of 20 |iM. The data represent the mean 
of net uptake rates of at least three independent experiments ± SE. The 
asterisks indicate significant differences between PLUTO mutants and 
the control based on Student's t-test (* = p<0,05; *** =p<0,005). 
doi:10.1 371/journal.pone.0091 343.g002 
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Table 2. Substrate affinity of PLUTO and PLUTO mutants. 







Mutant 


Uracil 


Guanine 


Adenine 




K„ ± SE [fiM] 


K„ ± SE [mM] 


K„ ± SE [nM] 


Control 


16.43 (±5.10) 


6.29 (±2.52) 


0.38 (±0.16) 


E227Q 






0.22 (±0.05) 


E227A 






0.21 (±0.05) 


E227D 


5.18 (±1.21) 


3.70 (±0.72) 


0.97 (±0.21) 


G147Q 


36.07 (±2,32) 


4.55 (±1,00) 


0.68 (±0.20) 


V145A G147Q 


20.15 (±3.73) 




0.81 (±0.17) 



Km" values are provided for uracil, guanine and adenine of different PLUTO 
mutants ± standard error (SE) and were determined after heterologous PLUTO 
expression in £ coti cells lacking the endogenous uracil importer uraA. 
Transport was corrected for uptake into uninduced cells. 
doi:l 0.1 371 /journal.pone.C091 343.t002 



Residues functioning as a selectivity filter 

In the PLUTO model, side chains of Ile-226, Phe-341 and Trp- 
223 are orientated upstream of the binding site, ideally positioned 
for defining the entrance of substrates in a trajectory leading to the 
binding site. Due to the strategic location of these residues above 
the binding site we hypothesize that these residues might be critical 
for determining selectivity for the physiological substrates of 
PLUTO by preventing substrate analogs from entering into the 
binding pocket. When the residues are substituted with short side 
chain residues they might allow passage of substrate analogues 
with bulky groups. Similarly, in the uric acid-xanthine permease of 
Aspergillus nidulans, a hydrophobic residue Phe-528 located above 
the substrate binding site has been proven to be a molecular filter 
[27]. To test our hypothesis, all three residues were mutated and 
subsequently tested for transport function. However, no differenc- 
es concerning direct uptake of radiolabelled uracil, guanine or 
adenine were observed (Table 1 , Figure 5), suggesting they do not 
play a role in substrate binding. 



In the second step, structurally similar substrates were tested in 
competition experiments. Adenine and guanine (20 |iM) transport 
was tested in the presence of 200 |tM hypoxanthine, allantoin, uric 
acid, benzylhydantoin, or thiamine. Uracil (20 |iM) transport was 
tested in the presence of 200 |tM thymine, thiamine, cytosine, 
orotate, dihydroorotate and dihydrouracil. Guanine transport in 
wild type PLUTO was affected by hypoxanthine and uric acid 
(Figure 5B) and uracil transport by thymine and orotate 
(Figure 5A). The latter two substances were also tested in direct 
uptake studies and were not identified as transport substrates. In 
W223A mutants, hypoxanthine and uric acid were no longer able 
to inhibit guanine transport (Figure 5B) while uracil transport was 
inhibited by thymine, thiamine and dihydrouracil (Fig 5A). In 
F341A and I226A mutants guanine transport was inhibited by 
hypoxanthine, no further effects of competing substrates were 
identified (Figure 5). No competitive inhibition of adenine 
transport was observed in any of the proteins (Figure 5C). 

Competitive binding of substrates 

Our PLUTO docking and site-directed mutagenesis results 
suggest the binding sites for uracil and guanine overlap whereas 
adenine adopts a different binding mode and binds deeper into the 
protein cavity. Therefore, we envisioned that competition studies 
with the substrates would provide some additional insight into the 
substrate binding. To check this, competition studies with all 
substrates were performed. Transport of 20 |tM ['*C]-uracil was 
reduced by adding 200 |iM guanine, but remained unaffected by 
the same concentration of adenine (Figure 6A). In case of guanine 
transport, both uracil and adenine had marked inhibitory effects 
(Figure 6B). Only adenine was not influenced by the presence of a 
ten times higher concentration of uracil or guanine (Figure 6C). 

When combining the obvious effects of competition studies and 
transport experiments with mutated PLUTO protein the following 
conclusions can be drawn: Uracil and adenine binding pockets are 
independent. Apart from the drastic mutations E227A, E227K 
and E227Q^ G147Q_ which affect all substrates, the remaining 
mutants affect either uracil and guanine or adenine and guanine 
transport (Figure 7). The same holds true for the competition 




Variable Conserved 



Figure 3. Position of functionally relevant residues in the PLUTO structural model. (A) The PLUTO homology model was built using I- 
TASSER server based on the open conformation of MHPI (PDB ID: 2JLN) and its quality was assessed using the PROCHECK program. Important 
transmembrane segments involved in nucleobase binding are marked in blue (TMs 1 a and 1 b), green (TM 3), pink (TM 8) and salmon (TMs 6a and 6b). 
Residues directly involved in substrate binding are VMS, G147, E227 and T425 (marked with yellow boxes and shown in the middle). (B) Evolutionary 
conservation profile of the PLUTO core protein (TM 1 -TM 1 0) calculated by Consurf server [23] shows a high conservation of the substrate binding 
pocket. The PLUTO homology model is colored according to the conservation score (based on ClustalW matrix). 
doi:1 0.1 371/journal.pone.0091 343.g003 
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A ^ ^ B 




uracil adenine 



Figure 4. Docking poses of substrates. (A) Schematic overview of the docl<ing poses of uracil (green), guanine (red) and adenine (blue). Uracil 
and guanine exhibit a similar binding mode, whereas adenine binds deep into the catalytic pocket of PLUTO. (B) Docking pose of guanine. (C) 
Docking pose of uracil. (D) Docking pose of adenine. Glide 5.0 (Schrodinger software) was used for docking purposes and the standard precision 
mode was selected to dock the substrates into the receptor site. 
doi:10.1371/journal.pone.0091343.g004 



experiments. Adenine transport is only afTected by the two double 
mutations: E227CiT425A and E227Q, G147Q, (Figure 7). Uracil 
transport is afiected by the mutations E227Q, V145A G147Q, and 
in addition by guanine as competitor (Figure 7). Guanine transport 
is affected by all mutations and in addition by both other substrates 
as competitors (Figure 7). 

Discussion 

A systematic mutational analysis of residues forming the 
putative substrate binding pocket in PLUTO has been performed 
to delineate their role as structural determinants for substrate 
binding. In the first step, conserved amino acid residues were 
identified using a multiple sequence alignment of NCSl proteins, 
including PLUTO and MHPl [4]. Two out of six residues 
identified as putative substrate binding site in the MHPl crystal 
structure differ compared to PLUTO. Besides differences in key 
amino acid residues, the substrates transported by PLUTO and 
MHPl are different. PLUTO transports uracil, adenine and 
guanine, whereas MHPl facilitates transport of 5-benzyl hydan- 
toin and 5-indolyl-methyl hydantoin [4,28]. When Glu-227 and 
Gly- 1 47 were exchanged to those amino acids present in MHP 1 in 
the double mutant E227Q,G147Q, PLUTO substrates were no 
longer transported (Figure 2 A-C). Importantly, we confu-med the 
presence of similar amounts of protein in the double mutant and 



control cells by western blot on total membrane preparations from 
induced E. coli cells (Figure SI). To test whether PLUTO can 
transport the known substrate of MHPl benzyl hydantoin uptake 
by wild-type PLUTO was studied using competitive binding 
analysis. However, no indication for a transport of benzyl 
hydantoin was obtained. These results indicate that although 
PLUTO and MHPl belong to same protein family, their active 
sites are arranged differently, thus defining a molecular determi- 
nant for substrate selectivity. 

To gain further insight into substrate binding, in a second step 
the homology model was used to dock all the three physiological 
substrates of PLUTO into the protein cavity to develop a three- 
dimensional model of PLUTO-substrate binding. Overall the 
homology model is similar to the crystal structure of the outward 
facing conformation of MHPl. Superimposition of all Co; atoms in 
PLUTO and MHPl gives a RMSD of 2.7 A. 4 TM helices 
including TMl, TM3, TM6 and TM8 form the ligand-binding 
site, in which TMl and TM6 show unwound regions, whereas 
TM8 is slightiy bent. A common arrangement with unwound and 
bent hehces in the hgand binding site is observed in both, MHPl 
and FCYB models. These structural conformations are important 
for transporter function, although one notable difference was 
found in the unwound region of TM6, which is more extended in 
PLUTO than in MHPl. We also found that some of the residues 
contributing to the structural organization at the ligand-binding 
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Figure 5. Competition studies with PLUTO mutants putatively acting as specificity filter. Direct uptake studies were performed after 
heterologous expression of PLUTO (black) and PLUTO mutants W223A (white), F341A (spotted) and I226A (striped) in f. coli cells lacking the 
endogenous uracil transporter uraA. The uptake was measured for uracil (A), guanine (B) and adenine (C) with concentrations of 20 [iM and several 
competitors were added with ten-fold excess (200 nlVI). The data represent the mean of net uptake rates of at least three independent experiments ± 
SE. The asterisks indicate significant differences between PLUTO mutants and the control based on Student's t-test (* =p<0,05; ** =p<0,01; *** 
= p<0,005). 

doi:l 0.1 371 /journal.pone.0091 343.g005 



site are not conserved between NCSl members, such structural 
heterogeneity might be a prerequisite for the substrate specificity 
within the family members. 

This information allowed us to identify key residues in ligand 
binding and to perform a comprehensive mutational mapping of 
these residues. Several important features regarding the molecular 
determinants of PLUTO substrate binding emerge from the study. 
First, the mutational analysis clearly identified Glu-227 as a critical 
determinant for guanine and uracil binding. Second, docking 
analysis identified Thr-425 in addition to Gly-147 as possible 
adenine interaction partner. However, although the binding of 
adenine is independent of Glu-227, the mutation E227Q_G147Q^ 
completely abolished the transport of adenine. Third, as evidenced 
by the substrate competition studies and docking pose, the binding 
sites for uracil and guanine overlap, whereas adenine binds deeper 
into the cavity. 

Some of the active residues showed backbone H-bonding with 
the substrates, however validation of such interactions is difiicult 
by single mutations. Therefore, double mutations of these residues 
with E227Q,were generated. E227Q_was chosen because here the 
differences in response of uracil, guanine and adenine were most 
pronounced. In fact, changing Thr-425 to alanine in E227Q_ 
T425A led to highly reduced adenine and guanine transport. This 
fits nicc-ly to the docking models of both purine bases, which 
interact with Thr-425. T425A alone showed no effect on adenine 
transport. This indicates that mutation of Glu-227 somehow 
modifies the binding pocket but the second mutation is required to 
restrict purine transport. V145A G147Q, most markedly affects 
guanine transport. Both residues bind guanine according to the 
docking studies, whereas adenine is only interacting with Val-145 
and uracil does not interact with any of the two. However, adenine 
transport is mainly influenced by modification of Gly-147 in close 
proximity to Val-145. Possibly our docking model does not allow 
exact prediction of H bond formation between adenine N7 and 
the corresponding amino acid residue. 

To gain further insight into the role of residue Glu-227, the 
mutants E227D, E227Q, E227K and E227A were constructed 
and analyzed. Reduction of the acidic chain of Glu-227 by one 
methylene (E227D) reduced transport of uracil and guanine. At 
the same time the affinity towards uracil increased in this mutant 
(Table 2). Interestingly, a similar mutation in the uric acid/ 
xanthine transporter UapA of A.sperff.llus nidulans E356D led to 
dramatically increased binding of the substrate but reduced 
transport [29]. Glu-227 is located in the middle of TM3, thereby 
residing between substrate binding-site and co-ion binding-site and 
possibly connecting these two sites. Strategic location of Glu-227 
might allow the residue to play a dual role in both substrate and 
ion binding, thus becoming a key element in coupling. A dual role 
of Glu-269 in both substrate and ion binding has already been 
observed in LacY permease [30,31]. Although LacY is only 
distantiy related to the bacterial nucleoside transporter NupG, 
both share common mechanistic features, as recentiy shown [32] . 
Moreover, Glu-264 aligning close to Glu-269 in LacY was shown 
to be important for uridine transport in NupG [32]. Probably due 
to its dual role in PLUTO, Glu-227 coordinates between ion and 
substrate binding. Once the residue is substituted with aspartate 



such link is broken, which might explain why E227D in spite of 
having higher substrate affinity showed reduced transport. 
Adenine transport is difierent, as the affinity of PLUTO is much 
higher with a Km value of 0.38 |J,M. Although the affinity drops in 
E227D to 0.97 |iM this is still 7-10 times higher compared to 
uracil or guanine. The high affmity of PLUTO for adenine is 
probably one reason for the low effect of mutations and 
competitors on adenine transport (Figure 5C, 6). In general, 
Glu-227 and Gly-147 mutants showed no obvious alterations in 
their substrate affmities. 

Recently, the purine/cytosine transporter FCYB from A. 
nidulans was analyzed by modeling, docking and mutational 
studies [8]. Like PLUTO and MHPl, FCYB is a member of tiie 
NCSl family. In FCYB five amino acid residues were identified as 
critical for substrate binding, these are Ser-85, Trp-159, Asn-163, 
Trp-259 and Asn-354. These align to Gly-147, Trp-223, Glu-227, 
Trp-342 and Asn-43() in the PLUTO sccjucnce, respectively 
(Figure lA). All of these PLUTO residues were mutated in the 
course of this work. One can differentiate between mutants 
directly involved in substrate binding (Ser-85, Asn-163, and Asn- 
354) and residues with a function in stabilizing substrate binding 
by pi-stacking interactions (Trp-159 and Trp-259) [8]. With 
respect to direct binding, only Glu-227 single mutants markedly 
affected the PLUTO transport capacity, equivalent to Asn-163 in 
FCYB. Single mutants of the other residues Gly-147 and Asn-430 
showed only weak or no effects (Table 1). The tryptophan residues 
involv(-d in /;!-stacking interaction in FCYB are also conserved in 
PLUTO. However, due to an antiparallel orientation in PLUTO, 
different to the situation in FCYB, the /)i-stacking interactions are 
weak and no effect on PLUTO facilitated transport of any 
substrate was identified in corresponding mutants. Thus, these 
interactions seem to be of minor importance for PLUTO function. 

It was concluded that residues corresponding to Trp-159, Trp- 
259 Asn-350, Pro-353 and Asn-354 in FCYB are major elements 
of the substrate binding site in all NCSl members [8]. As discussed 
above this does not hold true for both tryptophan residues in 
PLUTO. Asn-350 equivalent to Asn-426 in PLUTO was also 
mutated in the course of this work and shown to have no effect on 
transport. Furthermore, two additional residues Val- 1 45 and Thr- 
425 were identified as important components of the substrate 
binding site. These results highlight the difference between the 
substrate binding properties of PLUTO and FCYB. Moreover, 
Asn-430 in PLUTO is highly conserved in NCSl members and 
mutations extinguish the transport capacity of FCYB, but not of 
PLUTO [4,8]. 

Asn-163 is conserved in nearly all purine specific FCY-like 
NCSl proteins and is replaced by Glu-227 in PLUTO [4,8]. 
When this residue was mutated to Leu and Gin the transport 
capacity of FCYB was highly reduced or lost. Similarly, replacing 
Glu-227 by Gin exerted the strongest observed reduction in 
PLUTO mediated transport. Only Glu-227 seems to allow for a 
high capacity of purine and pyrimidine nucleobase transport of 
PLUTO. Interestingly, the yeast uracil transporter FUR4 carries a 
Gin residue at this position. Obviously, the important function of 
residue 227 awaits further clarification. 
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Figure 6. Substrate competition studies. Direct uptal<e studies 
were performed after heterologous expression of PLUTO in £ coli cells 
lacking the endogenous uracil transporter uraA. The uptake of uracil (A), 
guanine (B) and adenine (C) was measured after the incubation with 
20 nM substrate for 2 minutes (control). To investigate the effect of the 
other substrates, they were added to the uptake medium with ten-fold 
excess (200 nM). The data represent the mean of net uptake rates of at 
least three independent experiments ± SE. The asterisks indicate 
significant differences between PLUTO mutants and the control based 
on Student's t-test (* = p<0,05; *** = p<0,005). 
doi:1 0.1 371/journal.pone.0091 343.g006 

All FCYB substrates, adenine, hypoxanthine, guanine and 
cytosine interact with the same residues and a simple model of 
competition between substrates and the binding site is assumed 
[33]. This situation is obviously difTerent in PLUTO where the 



competition 



transport 




V145A Q) I© 



G147Q 

Figure 7. Schematic overview of competition between tKie 
different substrates of PLUTO (red) and transport activities of 
different PLUTO mutants (black). The presence of uracil leads to a 
competitive inhibition of guanine transport and vice versa. Further- 
more, adenine leads to a competitive inhibition of guanine transport. 
The mutation E227Q leads to a loss of uracil and guanine transport 
actitivity (-), whereas the double mutation V145A G147Q only affects 
guanine transport. In both double mutants E227Q T425A and G147Q 
E227Q guanine and adenine transport is affected. 
doi:1 0.1 371/journal.pone.0091 343.g007 

substrates interact with different residues and thus show different 
responses in competition experiments. 

A further aim was to identify the molecular basis of substrate 
selectivity using competition studies of functional mutants towards 
a plethora of substrates. In the PLUTO model Ile226, Phe341 and 
Trp223 are located above the binding site, which led to 
speculations about further roles of these residues in substrate 
selection. Substitution of the above mentioned residues with 
alanine does not affect transport and the specificity profile of 
PLUTO for its physiological substrates, suggesting that these 
residues are not an essential part of the substrate binding site. 
Among these mutants, W223A showed most dramatic results, as 
competition inhibition profiles of both uracil and guanine for non- 
physiological substrates were altered in the W223A mutant. Thus, 
Trp223 might be involved in the 'selection' of physiological 
substrates by performing pi-pi interactions with particular 
substrates and not with others, depending on how substrates are 
orientated in the nearby binding site. Furthermore, TrplI? in 
MHPl, equivalent to Trp223 in PLUTO, showed significant side 
chain movement during transition between outward to substrate- 
bound occluded conformation, indicating its role in either 
substrate binding or selection [18]. However the competition 
studies with mutant W223A do not provide strong evidence for a 
role in substrate-selectivity but might be useful for further 
investigations about the molecular basis of substrate selection in 
PLUTO and NCSl members in general. 

From the results obtained in uptake and competition studies as 
well as docking analysis, a model was generated (Figure 7). From 
this model, it can be concluded that guanine transport is most 
strongly affected by mutations and competes with uracil and 
adenine. This reflects the fact that guanine shares common amino 
acid residues with uracil (E227) and with adenine (T425 and 
G147) binding pockets. Furthermore, inhibition of guanine 
binding by excess of adenine can be due to the relatively large 
size of both substrates which could lead to steric clashes during the 
binding and eventually reduce the transport. Uracil only interacts 
with guanine as it shares no interacting amino acid residues with 
adenine. Finally, adenine is unaffected by competing substrates as 
it is located at a different position within the substrate binding 
pocket and exhibits the highest substrate affinity. 
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From the competition studies with purine related compounds, 
hypoxanthine was identified as a putative further substrate as it 
shows marked inhibition of guanine transport. Whereas the 
physiological function of uracil import into plastids is quite clear, it 
is unknown whether adenine, guanine or hypoxanthine might be 
further processed in this compartment. As discussed elsewhere, it is 
unlikely that adenine recycling to AMP by adenine phosphor- 
ibosyl-transferase (APRT) takes place in plastids [6]. Remarkably, 
mutants impaired in adenine recycling due to a mutation in the 
main APRT isoform showed a 50% increase in intracellular 
adenine levels accompanied by improved growth and stress 
tolerance [34]. The observed stress tolerance was attributed to a 
signaling function of adenine. However, it was not clarified in 
which intracellular compartment such signaling was exerted. It 
was shown that guanine is recycled to GMP by hypoxanthine/ 
guanine phosphoribosyltransferase which represents a single copy 
gene in Arabidopsis. Corresponding knockout mutants were 
almost completely resistant to toxic 8-azaguanine during germi- 
nation which is in line with the inability of such mutants to 
incorporate 8-azaguanine into RNA and DNA, leading to its 
toxicity [35]. There is no indication that this enzyme or 
Arabidopsis APRTs are located to plastids by bioinformatic 
analysis [6] or analysis of corresponding mutants [35] nor is the 
catabolism of purine nucleobases located in plastids [36] . 

Conclusion 

In conclusion, it was shown that by combining structural and 
functional studies including mutants, critical amino acid residues 
involved in interaction with the different substrates of PLUTO 
could be identified. By this, a homology model constructed on the 
basis of the structural information from MHPl was validated. 
Significant differences in structure function relations between 
PLUTO and the only other characterized eukaryotic member of 
the NCSl family, FCYB were identified. Although the PLUTO 
structural model predicts a single binding site, not all substrates 
compete with each other due to minor differences in the 
interaction with specific amino acid residues. This illustrates that 
competition studies often cannot replace direct transport studies. 
Another prominent example for this is the sucrose transporter 
SUC5 facilitating sucrose and biotin transport without reciprocal 
inhibition [37]. These results highlight key differences between 
active site architecture of PLUTO and NCSl family members, 
which will be important for the prediction of the substrate 
specificity in other NCSl proteins. 
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